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' A theoretical analysis of the ground state of long-chain (K@C6o)n is presented. Within mean field 

. theory, a ferroelectric ground state is found to be stable because of the pseudo-Jahn- Teller mixing 

^ ■ of the b\u and the 629 band with a zone-center optical phonon involving the displacement of the 

^ ' endohedral ions. A phase diagram for this model is derived in the narrow bandwidth regime. 
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■ Introduction 

• 1— I \ 

^ • FuUerenes have proven to be useful building blocks in the synthesis of novel molecular solids. Metals, insulators, 
' I , ferromagnets, and superconductors have all been fashioned out of fuUerenes. A variety of fuUerene-based polymeric 
• structures have also been synthesized and characterized |[ ^ ^, |^. Linear chains of CgQ have been produced jl) by 
' the coevaporation of K with Cgo- The external K atoms readily donate their valence electrons to the Cgo molecules, 

forming a charged "pearl necklace" of Ceo anions (see Fig. |l|). 
' Endohedral doping, where one or more atoms are positioned inside the Cgo cage, offers an alternative way to modify 

a the properties of fuUeride solids. Theoretical studies show that sizeable electric dipole moments]^, ^, ^, |l^, 
I ' can result from the endohedral doping of isolated fuUerenes with metal atoms, so-called metallofuUerenes. The dipole 
' formation can be understood pl]| as a consequence of the pseudo-Jahn- Teller effect (PJT) where the ti„ and tig 
^ ' electronic states of the metallofuUerene are mixed by a tiu distortion involving the displacement of the central metal 
O \ atom, the {tiu (Btig)® ii„ model. 

The properties of a linear chain of metallofuUerenes are investigated in what follows. A lattice generalization of the 
{hiu ® b2g) ® bsu PJT model is explored. A phase diagram for this model is derived in the narrow bandwidth regime. 
It is shown that within mean field theory, a ferroelectric ground state is stable when the coupling constant g exceeds 
^ ' a critical value go which depends on the cation force constant k and the biu — 629 energy spacing A. 
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Model Hamiltonian 



' ESR measurements on exohedral polymeric (KC6o)n are consistent with one electron transfer from K atoms to the 
Cgo chains. Adjacent Cgo molecules bond with a preferred orientation as indicated in Fig. |^. The axis of the chain 

■ ' bisects a pair of hexagon-hexagon bonds on every Ceo molecule. It is assumed that endohedral polymeric (K@Ceo)n 
will have the same structure. 

Such an arrangement of metallofuUerenes with the ions centered with respect to the Ceo cages will have a point 
group symmetry of T)2h, as the chain axis is coincident with one of the molecule's two- fold rotation axes. Facing 
rH • hexagon- hexagon bonds on adjacent Ceo molecules align, and the two intermolecular bonds together with the two 
Q facing hexagon-hexagon bonds complete a four-membered ring. A coordinate system is chosen so that the chain axis 
O is labeled by z and the four-membered rings lie in the yz-plane as indicated in Fig. |l|. 

^ _ In an isolated alkali-doped metallofuUerene, the highest occupied molecular orbital (HOMO) has a three-fold spatial 
degeneracy and tranforms as ti„ under the icosahedral group. With the site symmetry lowered to D2h, this three- fold 
degeneracy is lifted and the ti„ multiplet is split into its components under D2h- ^lu, ^2u, and h^u- Ligand-field 
. analysis indicates that the h2u and b3„ levels which transform as y and x will be closely spaced in energy, while the 

■ - ' hiu will be lowered by more than the splitting S between b2« and b3„. 

It can be argued using symmetry considerations that the h2u state will be lower in energy than the b3„ state. 
The HOMO is composed primarily of radially directed atomic p orbitals. Since the bsu state must transform as x, 
the yz-plane must be a nodal plane for this state, as any combination of pz orbitals centered on the sites of the 
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four-membered ring would be even under reflection through the j/z-plane. Thus the C atoms participating in the 
hexagon-hexagon bonds on axis can not contribute Pz character to the ba^ orbital. 

In the b2u state, however, these C atoms can contribute Pz character and satisfy the symmetry requirement of 
oddness under y — > —y. Consequently, the state has a lower energy through cr-bonding to nearest neighbors along 
the chain axis. 

The lowest unoccupied molecular orbital (LUMO) in an alkali-doped metallofullerene transforms as tig. In analogy 
with the HOMO, in the polymer, this niultiplet is split into 623, and &3g with the ordering given in Fig. ^. 

It has been previously shown |pl| that there is strong vibronic mixing of the tiu and tig orbitals in metallofuUerenes 
through central atom displacements. A lattice generalization of that pseudo-Jahn- Teller (PJT) model is now described. 

The Hamiltonian is given by 

H — Ho + TYpjT + Holast (1) 

where 

^0 = -t^{c\„C^n+l + 4„C2„+1 + ll-C-) + ^4„C2„ (2) 
n n 

where n runs over endofuUerene site, cj„ (c^^) creates (annihilates) an electron in orbital biu at site n, (c2„) 
creates (annihilates) an electron in orbital b2g at site n, and A is the energy spacing between the biu and &2g orbitals. 
The hopping matrix is taken to be diagonal in orbital type and spin, and the spin indices are suppressed for clarity. 
The pseudo-Jahn- Teller Hamiltonian on the ?ith endofuUerene is taken to be 

HpjT = - — 5I(cI«C2„Q« + h.c.) (3) 
a ^ — ' 

n 

where a is the equilibrium spacing between successive endofuUerenes, and g is a dimensionless coupling constant. 

The point group symmetry only allows for vibronic mixing of biu and 629 orbitals with a displacement transforming 
as &3„. Hence, Qn is a displacement on the nth endofuUerene in the x direction, perpendicular to the plane of the 
four-membered rings. 

'^eiast is the elastic energy for the displacements of the K+ ions. 

"Hclast = 2 X! (4) 




FIG. 1: Bonding structure of (KCeo) 
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FIG. 2: Correlation diagram of tin and tig manifolds under an axial ligand field. 



With the following transformations, 



(5) 



the Hamiltonian can be written in the fc-space representation 



(6) 



(7) 



where eg{k) = A — 2tcosk and e„(fc) = — 2tcosfc. 

In analogy with (KC6o)n, the ti„ and ti„ bands in (K@C6o)n are expected to have a width of the order of i eV. 
Hence, (K@C6o)n is expected to be in the narrow bandwidth regime where the 6i„ — 62^ orbital splitting is much 
larger than the bandwidth, A » it. 
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FIG. 3: Ground state energy per unit length E vs. Q = 2gA\Qq\/a\/ N for q = and q = tv (dashed). The values g = 20, 
t = 0.5 eV, A = 1 eV and k, = 1.6 eV/A^ were chosen. Top right graph shows detail in vicinity of a minimum. 



Mean Field Theory 



With a single mode of distortion given by a non- vanishing Qq, TL can be diagonalized with a canonical transformation 
so that 



(8) 



where 



Ek± = ^(e«(fc) + es(fc + q)) ± + q) - ^uiW + 



(9) 



and nk± are the usual occupation number operators of the vibronic quasiparticles. 

The ground state energy of chain with one electron per endofullerene on average is obtained by half-filling the lowest 
energy vibronic band. In the continuum limit, 



J-7v/2a ^ 



(10) 



The ground state energy can be calculated to C'(-^)^ 



NaA 



l-(f)fcos2W2-v/l + Q?-(f)f^ 



' 4t\2Ql sin^ qa/2 



4(1+Q?)2 



(1 + i sinqa) 



+ «|Q,|V2 



(11) 



where Ql = Ag^\Qq\^ /Na"^. 

At g = 0, the energy is a minimum for non- vanishing Qg when g > go where go = Ka^/2A. The results for the 
ferroelectric and antiferroelectric cases are displayed in Fig. ^. 
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FIG. 4: Zero-temperature phase diagram for (K@C6o)n in the narrow bandwidth regime u = ^ << 1. 

When g > go the energy minimum is found at g = 0, corresponding to an ferrodistortive transition. This displace- 
ment creates local electric dipoles on the endofullerenes that are aligned perpendicular to the yz-plane, yielding a 
uniaxial ferroelectric phase (FE). A paraelectric (PE) phase exists along the line segment g > go and t — where 
independent, local dipole moments form on the endofullerenes. The symmetric phase (S), where no local dipole 
moments form, is stable for g < go. These results are summarized in zero-temperature phase diagram pictured in 
Fig. |. 
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